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Abstract. 

In quantum information theory, the reliable and effective detection of entanglement 
is of paramount importance. However, given an unknown state, assessing its 
entanglement is a challenging task. To attack this problem, we investigate the use of 
random local measurements, from which entanglement witnesses are then constructed 
via semidefinite programming methods. We propose a scheme of successively increasing 
the number of measurements until the presence of entanglement can be unambiguously 
concluded, and investigate its performance in various examples. 
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1. Introduction 

Entanglement is a valuable resource for accomplishing quantum information tasks 
with a performance exceeding what can be optimally achieved classically. Thus, 
the unambiguous detection of entanglement is a fundamental problem in quantum 
information processing. However, even though in the limit of high Hilbert space 
dimension, a randomly chosen quantum state is almost certainly close to maximally 
entangled, this is in general a hard task [1], 

Nevertheless, important insights have been gained for the case that the state to be 
analysed is at least partially known. Given some prior knowledge of the density matrix 
(obtained, for example, by assuming that the given state is close to a target state that 
has been imperfectly prepared), few local measurements suffice to detect entanglement 
[2], In a different vein, local measurements along randomly chosen directions may be 
used to certify nonclassical correlations by means of a Bell-inequality violation, if the 
state distributed between the different parties is known, thus obviating the need for a 
shared reference frame [3, 4, 5, 6]. 

How can entanglement be detected in the case in which one does not have any 
knowledge of the given state? A method to characterize entanglement in such a case 
is given by quantum tomography [7], i.e., reconstructing the full density matrix from a 
pre-chosen set of measurements. 

However, this method has several drawbacks: first of all, it is enormously resource¬ 
intensive, with the number of measurements that need to be performed scaling 
exponentially with the number of particles of the system. Furthermore, recent results 
have shown that the most commonly used tomography schemes may be subject to 
systematic errors that can lead to overestimating the entanglement contained in a 
quantum state [8]. 

Additionally, if one is merely interested in the question of whether a given state is 
entangled, full tomography yields an excess of information; thus, it may be feasible to 
extract the required knowledge using fewer measurements. An approach in this direction 
was proposed in Ref. [9]. 

In a similar vein, in Ref. [10], a procedure for measuring the entanglement of an 
unknown state by successively measuring witness operators was proposed, which was 
experimentally realized in Ref. [11]. However, the approach in Ref. [10] is not very 
effective for mixed states: as shown in Ref. [11], about 67% of mixed states still require 
a tomographically complete set of measurements. An additional adaptive scheme is 
presented in Ref. [11], which considerably improves the performance for mixed states; 
however, only the two-qubit case is considered there. 

Furthermore, in contrast to these works, in this article, we investigate the possibility 
of detecting the entanglement of an unknown state using random measurements. Our 
conclusions thus do not necessitate the performance of any fixed measurement strategy. 
The intuition behind this is that if one is given an unknown, i.e. random, state, 
there is no canonical way to choose measurements tailored to optimally determining 


Detecting Entanglement of Unknown Quantum States with Random, Measurements 3 

the entanglement of this state. Indeed, any particular fixed choice of measurement 
sequence will introduce a bias, leading to certain states being systematically detected 
earlier than others, which is avoided by a random choice. 

The method of choice for entanglement detection we will be concerned with is that 
of so-called entanglement witnesses [12, 13, 14]. An entanglement witness is a Hermitian 
operator W such that 

tr (Wa) > 0 ( 1 ) 

for all separable states a, and 

tr (Wp) < 0 ( 2 ) 

for at least one entangled state p. Thus, a negative expectation value as a result of 
measuring W certifies the presence of entanglement. 

Note that the absolute of a negative expectation value for W also leads to a lower 
bound on an entanglement measure [15]. 

For simplicity, we will in the following mainly concentrate on a certain class of 
entanglement witnesses, the so-called decomposable witnesses [12, 16]. These are witness 
operators W that can be written in the form 

W = Q Ta + P. (3) 

Here, the superscript T 4 denotes the partial transposition operation with respect to 
subsystem A, and Q and P are both positive semidefinite operators, that is, they 
have no negative eigenvalues. This is the most general witness capable of detecting 
states with nonpositive partial transpose: tr (pW) < 0 implies p TA ^ 0, since 
tr(pPF) = tr(Pp) +tr(Qp TA ), which can only be smaller than zero if p TA is not positive. 
Nonpositivity of the partial transpose is a sufficient criterion for entanglement in all 
dimensions, and is even necessary in case dim(p) <2x3 [12], which yields the famous 
PPT-criterion for entanglement detection. 

The paper is organized as follows. In Section 2, we investigate the possibility of 
detecting entanglement using a randomly chosen witness operator. Then, in Section 3, 
we focus on constructing witnesses from random local measurements. To this end, in 
Section 3.1 we employ a decomposition of witness operators, and carry out optimization 
by means of a semidefinite program in order to find the best witness constructible from 
a given set of measurement settings and outcomes. Then, in Section 3.2, we propose a 
scheme that successively increases the number of measurements, until the semidefinite 
program certifies the presence of entanglement. In Section 4, we test this scheme on 
various experimentally relevant examples. Section 5 is devoted to investigating the 
statistical robustness of our method, and in Section 6 , we conclude. 

2. Randomly Generated Entanglement Witnesses 

As a first attempt towards addressing the problem of unknown-state entanglement 
detection, one might, reasoning that since there is no way to single out a particular 
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witness, simply consider using witnesses that themselves are (Haar-)randomly drawn 
from the set of all Hermitian operators acting on RLab = 'hiA Operationally, this 

corresponds to carrying out a random global measurement on the system. 

The first step, then, is to identify witness operators among all possible 
measurements. This can be done by checking two criteria which both have to be fulfilled: 
First, one checks whether the operator is indefinite. If it is not, then it cannot possibly be 
a witness, as its expectation value would have the same sign on separable and entangled 
states alike. Second, one checks whether the operator’s partial transpose is positive: 
Since for any separable state cr, o Ta is again a valid (separable) density operator, we 
can find o' such that o' Ta = o, and thus, if W Ta > 0, 


tr(HTx) = tr(Wo ,TA ) = tr (W Ta o') > 0, 


( 4 ) 


meaning that the expectation value of W is positive on all separable states. Hence, as 
W is indefinite by the previous requirement, it must detect at least one entangled state. 

In general, not all entanglement witnesses can be found in this way. The reason 
for this is the existence of entangled states in Hilbert spaces of dimension greater than 
d = 2x3, whose partial transposition is positive, but which are nevertheless entangled, 
and hence, cannot be detected by the above method. For the remainder, an optimization 
procedure that minimizes the expectation value of the operator on the set of separable 
states, based on the overlap minimization algorithm presented in Ref. [17], will be 
used. Positivity of this minimal expectation value then certifies the operator as an 
entanglement witness, provided a global minimum is attained. To certify that this 
occurs with a high degree of confidence, the optimization is performed 10 3 times with 
different initial conditions. 

In order to test the witnesses obtained in this way, we generate random states, and 
attempt to detect them, using randomly chosen operators. Random (mixed) states are 
produced as follows. First, a random unitary matrix can be generated by orthogonalizing 
a matrix with uniformly random, bounded complex entries using the Gram-Schmidt 
method (which guarantees Haar-randonmess [18]). A random k x n —dimensional pure 
state is then obtained by generating a random unitary and extracting one of its column 
vectors. Then, the k —dimensional environment is traced out, yielding an n —dimensional 
mixed state randomly distributed according to [19] 


dji(p) oc @(p)5(tr(p) — l)detp fc n , 


( 5 ) 


where the step function 0 enforces positivity, and the Dirac 5 guarantees the 
normalization of the density matrix thus obtained. 

The use of this prescription is motivated by the fact that mixed states, in 
quantum theory, always arise from interactions with an environment, together with 
the assumption that the total (pure) system-environment state is distributed according 
to the unitarily invariant (Haar-) measure [19]. 

Additionally, random maximally entangled states can be generated by starting from 
a given d-dimensional maximally entangled state, | v b + ) = ^ JA | i) 0d , and transforming 
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it with random local unitaries Ui [19], i.e. 


|Vv) = (g)^|* + }, 

i 


( 6 ) 


where the Ui are random unitaries acting on the local Hilbert spaces 'H ll thus effecting 
a random transformation leaving the amount of entanglement constant. 

However, the results of carrying out this investigation are discouraging: Generating 
random Hermitian operators by generating a diagonal matrix D with uniformly random 
real entries sampled from a bounded interval, and then forming 


H 


UDW 
tr(UDUQ ’ 


(7) 


where U is a random unitary matrix, we find that even for a two-qubit Hilbert space, 
witnesses are rare among randomly chosen operators. Indeed, among 10 5 candidate 
operators, only a fraction of 1.73±0.05% were identified as witnesses by the positivity of 
their partial transpose, with a further 1.34±0.04% being found to be witnesses by means 
of the overlap-minimization algorithm. Here and in the remainder of this section, the 
specified uncertainties are due to finite-size statistical effects. Even further discouraging 
is the fact that typically, witnesses being chosen according to this prescription have 
a very small likelihood of actually detecting the state one is presently interested in— 
testing witnesses found according to both criteria on a set of random states, it was found 
that only a fraction of (9.3±0.7)-10~ 6 of states was detected by witnesses having positive 
partial transpose, while a further (13 ± 1) • 10“ 6 were detected by witnesses found using 
the overlap-minimization approach. Even in the case of maximally entangled states, 
only (1.094 ±0.007) • 1CT 2 , respectively (1.092 ± 0.008) ■ 1CT 2 , of all states were detected 
using these methods. 

As a first approach towards increasing the efficiency of this method, one can 
consider increasing the number of witnesses generated. Since we do not want to perform 
additional measurements, the only way to do so is to add a (positive or negative) multiple 
of the identity onto a randomly generated observable that does not meet our criteria. 
While it is always possible to create an indefinite operator in this way, positivity of its 
partial transpose cannot always be guaranteed. Thus, whether such a shift leads to a 
witness depends on the eigenvalue spectrum of the original random operator and its 
partial transpose. 

Additionally, it is possible to optimize witness operators by means of adding a 
(negative) multiple of the identity, that is, obtain a new witness operator capable of 
detecting additional states. However, the gains of such a procedure are slight, and 
numerically, the witness detection efficiency remains of the same order of magnitude. 

Fortunately, as we will show in the remainder of the paper, this is not yet grounds 
to give up on the project of detecting entanglement using random measurements. In the 
following section we will introduce the method of constructing entanglement witnesses 
from random measurements carried out on the local Hilbert spaces. 
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3. Entanglement Witnesses from Random Local Measurements 

From this point on, we will only consider decomposable witnesses, that is, witnesses of 
the form given in Eq. (3). The reason for this is, as we will see, that such witnesses can 
be efficiently constructed using a semidehnite program. Hence, from now on, the term 
‘witness’ should always be understood to mean ‘decomposable witness’, unless specified 
otherwise. 


3.1. Constructing Witnesses 

Any Hermitian operator O in BQHab ), the space of linear bounded operators acting 
on the total Hilbert space Ear, can be decomposed into local Hermitian operators 
Ai G B(Ea) and Bj e B{Hb) such that 

o = CjjAj 0 Bj , ( 8 ) 

ij 

with Cij € M. Thus, Ending a witness operator W, given local measurements Ai and Bj, 
amounts to finding coefficients such that 

tr(Wp) = ^2 Cijtr(Ai 0 Bjp) < 0, (9) 

ij 


for some entangled p, and 


tr(W p) = '22 c ij^(A i 0 Bjp ) > 0, 
ij 


( 10 ) 


for all separable p. This translates to the following optimization problem: 


minimize: c • m 

subject to: W = J2ij c ijAi 0 Bj 
W = P + Q Ta 
P> 0 
Q >0 

tr(W) = 1 

Here, c is the vector of the coefficients c^, and m = ((A* 0 Bj)) is the vector of the 
expectation values of Ai 0 Bj. Both vectors are obtained by mapping the indices i and 
j to a single index a. The semidefinite program (SDP) thus formulated calculates 
the operator W, such that W has the lowest possible expectation value given the 
experimentally obtained data, while being a normalized decomposable entanglement 
witness. If this expectation value now is negative, it can be unambiguously concluded 
that p is entangled. Note that the requirement of unit trace here does not provide an 
additional constraint on the witnesses found, but merely ensures their normalization: 
observables with vanishing trace cannot be witnesses, since their expectation value with 
the identity necessarily vanishes. Thus, the associated hyperplane contains the identity; 
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however, there always exists a ball of separable states around the identity, and hence, 
any such operator necessarily possesses a negative expectation value on some separable 
states. 

This raises the question of how many measurements need to be performed in order 
to be able to certify the presence of entanglement—or, in other words, how many 
measurements one needs minimally to be able to construct an entanglement witness 
for an arbitrary state. 

An upper bound for this number is given by the number of measurements in 
a tomographically complete set: since randomly drawn measurement directions are 
linearly independent, such a set forms a basis for the Hilbert-Schmidt space of operators, 
and thus, any operator can be written as a linear combination of them. 

For a system of dimension d = 3 a • d B , the density matrix p contains d 2 — 1 
independent parameters, each of which needs to be fixed by a suitable measurement. 
However, in total, (d\ — 1) • (d 2 B — 1) measurements of the form A; <g> Bj suffice, since 
one can infer from their outcomes the d 2 A — 1 + d 2 B — 1 additional values corresponding 
to the observables A* 8 1 and 1 8 -Bp 

3.2. Improving Detection by Increasing the Number of Measurements 

An immediate advantage of the method proposed above is that it can be used to devise 
a scheme which always terminates with the successful detection of entanglement, if 
entanglement that can be detected using the PPT-criterion is indeed present. This 
can be achieved by successively adding measurement directions on Alice’s and Bob’s 
local systems, then running the semidefinite program, until entanglement has been 
detected. In the worst case, the algorithm finishes when a tomographically complete set 
of measurements has been performed. 

In order to better illustrate this strategy, we order the measurements as follows. 

Table 1. Table of local measurement operators on Ha ®Hb- 


1a ® 1 B 

1a® Bi 

1a ® B 2 

1a ® B :i ••• 

A \ 8 1b 

Ai ® Bi 

Ai ® B2 

Ai ® H3 

A2 8 1b 

A2 ® B\ 

A2 ® B‘2 

A2 ® -B3 ■ ■ ■ 

A3 0 1b 

A3 8 B\ 

A3 8) B‘2 

A 3 8 B 3 ■ ■ ■ 


In this ordering, the measurements in the top row and leftmost column need not 
be performed, but instead come ‘for free’ upon performance of the measurements in the 
second row and second column. 

Different ways of traversing these measurements now suggest themselves. Strategies 
may differ in the experimental complexity they necessitate: We clearly want to minimize 
the number of measurements that need to be carried out in order to detect a given state; 
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however, at the same time, resetting measurement apparata is itself an investment of 
time resources and may also introduce a source of errors. If, for instance, we choose 
to first implement all measurements of the form A\ ® Bi, that is, leave the apparatus 
at A fixed, and only change the setting at B, we have to re-set the apparatus at B in 
order to measure the observables A 2 ® B i: A 3 ® Bi, and so on; but this is in general only 
imperfectly possible. 

Among the strategies to traverse Table 1, we sketch three obvious possibilities: 

(i) At each step, we add a new measurement, alternating between Alice’s and Bob’s 
detectors; then, we measure all combinations Ai <E> Bj of the measurements in the 
set. This corresponds to successively traversing the complete square of observables 
in Table 1. 

(ii) At each step, we add a new measurement on both sides, performing a new global 
measurement Aj 0 Bi. This corresponds to traversing Table 1 along the main 
diagonal. 

(iii) At each step, alternatively either Alice or Bob changes their measurement, while the 
other party continues theirs. Thus, we get the measurement sequence Ai <g) B\ —> 
A\ ® B 2 — ^ A 2 $5 B 2 —>-.... 

An intrinsic advantage is provided by strategy 2: since every new round conveys 
a maximum of new information, we may expect it to perform, on average, better than 
the other strategies; in fact, since one needs at least two distinct local measurements to 
detect entanglement [20], this is the only strategy that enables entanglement detection 
already in the second round. Thus, in the reminder of this paper, we will use this 
strategy. 

4. Detecting NPT-Entanglement of Unknown States 

In the previous section, we have established a method to detect entanglement of states 
about which nothing is known except the dimension, and hence, which can be assumed 
to be drawn randomly. In order to now put this method to the test, we will simulate 
its performance on states chosen randomly as described above. 

To quantify entanglement, we resort to the negativity, a well-known entanglement 
measure that was introduced in Ref. [22]. This is also motivated by the fact that we 
detect states exactly if their partial transpose is nonpositive, that is, their negativity 
is nonzero. Furthermore, the expectation value of (optimal) decomposable witnesses is 
equal to the negativity [21], 

The negativity is defined as 



( 11 ) 


where ||X||i = tr(\/XXQ denotes the trace norm, and the c Oi are the eigenvalues of p TA . 
Thus, the negativity of p is the absolute sum of the negative eigenvalues of its partial 
transposition [22], 
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The distribution of negativities is in general such that it is very unlikely to randomly 
draw a highly entangled state, and, even though the bulk of states gets more entangled 
in higher dimensions, this gets worse with increased system size. 

Our analysis concerns the performance of our method in detecting quantum states 
uniformly randomly chosen as described above. As a figure of merit, we use the number 
of measurements needed until a given entangled state is detected; thus, our algorithm 
first computes the expectation values of a certain number of random measurements A* 
and Bj given the random state p, and then executes the semidehnite program to see 
whether the state’s entanglement has been detected. If that is not the case, additional 
measurements are simulated one at a time, where each additional measurement is 
composed from local measurements in the form A* <g) Bj, and chosen according to the 
scheme discussed in Section 3.2. The algorithm is run until entanglement is detected. 
The results, for a system of two qubits, are shown in Figure 1, while Figure 2 analogously 
shows the same analysis for two qutrits. 

As one would expect, the method becomes more effective for more strongly 
entangled states, since these lie closer to the border of the set of all quantum states, 
and thus, less information is needed to distinguish between them and separable states. 
Compared with the naive method discussed in Section 2, we also observe a marked 
increase in efficiency: where the mere random measuring of entanglement witnesses 
detected even highly entangled states with infeasibly low probability, here, highly 
entangled states are detected after performing only a few measurements. 

For 3 x 3-dimensional systems, the gain in efficiency with more highly entangled 
states is even more pronounced. Although few states of truly high entanglement content 
are produced by our random method, even moderately entangled states are typically 
detected after only around 20 measurements, whereas full tomography would need 64 
measurements in general. 

Indeed, as shown in Figure 3, for maximally entangled 3 x 3-dimensional states, 
one can typically detect their entanglement after performing only 10 ± 3 measurements 
(where the uncertainty is due to finite size statistics). 

5. Statistical Analysis 

Up to this point, we have concentrated on the ideal case, that is, used the exact 
expectation values (A, ® Bj) in order to construct witness operators. However, in a real 
experiment, we do not have access to these values, but only to experimental data subject 
to statistical uncertainties and systematic errors. In the following, we will concentrate 
exclusively on statistical fluctuations; for a discussion of systematic errors in the context 
of detecting entanglement using witness operators, and a method to ameliorate these, 
see [23]. These fluctuations introduce an uncertainty into the experimentally obtained 
expectation values. Given some dichotomic operator M , if we perform N experiments, 
obtaining n + times +1 and n_ times the value —1, we can calculate the approximate 


Detected Fraction 
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Figure 1. (color online) 10 5 runs of the algorithm with 2 x 2-dimensional states 
drawn randomly. The data is normalized with respect to the binning of the negativity 
values, such that for each value of the negativity, it sums to one. For separable states 
or states of little entanglement content, it can be seen that the number of necessary 
measurements reaches the tomographically complete maximum, 9. 


expectation value 

M=i( „ + -„-) = i(2„ + -l), (12) 

where the bar denotes the experimentally obtained mean value. 

The value n + now carries a statistical uncertainty A n + = \jNp( 1 — p ), where p 
is the probability of obtaining the outcome +1 in a single measurement, since it is 
binomially distributed. Using standard error propagation, the resulting uncertainty of 
M is 

Aff = ^ An+ = 7y v/p(1 ~ p) ' (13) 

Of course, without knowing the state p, we do not have access to the probability 
p. However, as p(l — p) assumes its maximum of | at p = |, we can use the worst-case 
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Figure 2. (color online) 3 • 10 4 runs of the algorithm with 3 x 3-dimensional 
states drawn randomly. The data is normalized with respect to the binning of the 
negativity values, such that for each value of the negativity, it sums to one. Note 
that, since low-negativity states are unlikely to be generated by random draw, states 
with too little entanglement content are not shown; otherwise, the maximum number 
of measurements, reached at least in the case of separable states, would be 64, as 
required for tomographic completeness. 


approximation 

AM < —L. (14) 

~ y/N 

According to the second strategy defined in Section 3, the witness is of the form 

w = Yl c * M " ( 15 ) 

i 

with Mi = Ai® Bi. Here, we have the fortunate case that we need not worry about the 
covariance between different observables, as all are chosen independently. It would now 
be tempting to simply proceed using, again, error propagation, obtaining the formula 


AIM = . 

N 

However, doing so would neglect the fact that the coefficients c,; are not independent of 
the mean values Mj, since they are, in fact, derived from them. 

This difficulty can be overcome by dividing the data into two sets, using one to 
derive, via the semidefinite program, the coefficients c t , and then evaluating this witness 
using the expectation values generated by the other set of data (cf. [24]). That way, the 
coefficients are independent of the uncertainties in the expectation values of the second 
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Figure 3. (color online) 2 • 10 4 runs of the algorithm with maximally entangled 3x3- 
dimensional states drawn randomly according to Eq. 6. 


set, and thus, can be assumed to be simple constants, and Eq. 16 can be used to finally 
obtain 


AW 



< 



(U) 


In case every observable is measured an equal number of times, 
this further simplifies to 


AW < 



i.e. 


Ni 


N for all i, 
(18) 


This estimate can be used to facilitate the decision if, in an experiment, one should 
rather consider adding new measurements, or performing additional repetitions of the 
measurements already made: Especially in the case of an observed expectation value 
close to 0, a very small uncertainty is necessary to conclude the presence of entanglement; 
however, since additional measurements will tend to drive the expectation value closer to 
the state’s negativity, adding another measurement instead may be more advantageous. 

Consider, to this end, the example in Figure 4: There, for a single detection of 
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a low-negativity (A f = 0.0163) random two-qubit state, the maximum of the 3<r-error 
interval is plotted as a function of the number N of measurement repetitions. As can 
be seen, while a detection, that is, a certification of tr (Wp) < 0 with 3cr-conhdence, 
is possible for 6 global measurements, doing so requires a high number of repetitions 
for each individual measurement; however, adding another measurement, this number 
drops drastically. 



Figure 4. (color online) The maximum of the 3<r-confidence interval for the witness 
expectation value as obtained by the statistical estimate according to Eq. 18. 


6. Conclusions 

We have proposed a new method to detect the entanglement of quantum states, about 
which nothing is known except the dimension, using random local measurements. From 
these measurements, a witness operator is constructed by means of a semidehnite 
optimization procedure, minimizing the expectation value in a way that ensures that 
the operator is a valid witness at all times. Thus, once this expectation value becomes 
negative, one can unambiguously conclude that the state must be entangled. 

Our main objective was to investigate the feasibility of this method in an 
experimental context, and in particular, compare it to existing methods of investigating 
unknown states, such as for instance local tomography. To this end, we have carried out 
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numerical simulations, investigating the performance of our method using states drawn 
uniformly at random. 

We found a significant reduction in the number of measurements one needs to 
perform in order to certify the presence of entanglement, as compared to doing full 
tomography, yet still without introducing any additional assumptions about the state. 
For instance, in the case of a 3 x 3-dimensional system, maximally entangled states may 
be detected after only 10±3 random measurements, as compared to the 64 measurements 
needed for full tomography. 

An additional question is the performance of the method in higher-dimensional 
cases in comparison to tomography. Indeed, one might worry that the complexity of 
the SDP scales itself so badly as to effectively spoil any advantage gained due to the 
smaller number of measurements that need to be performed. 

However, while it is difficult to gather sufficient statistics for high-dimensional 
systems, the scaling of our method is encouraging: Even on a standard desktop system, 
a single instance of the optimization can be performed in a few minutes for an 8x8- 
dimensional system, while the performance of the required number measurements for 
full tomography is certainly on the edge of feasibility. Going to yet larger systems, in 
the case of d = 11 x 11, the optimization procedure takes about two hours, again on a 
simple desktop system. 

A question of special relevance to the experimental implementation of this method is 
the statistics needed in order to certify entanglement to some given degree of confidence. 
Here, too, our method has advantages compared to tomographic schemes, since we are 
able to quantify the degree of confidence with which entanglement has been detected, 
while many current tomography protocols suffer from the possibility of false positives 
([8]). In order to facilitate this, we derived an upper bound on the uncertainty of the 
witness expectation value, and provided a method to gauge the statistical fluctuations 
during an experiment, by binning the data into independent sets, and cross-checking 
the witness constructed from one set using data from the other [24], Simulation of this 
procedure showed a good robustness of the method to statistical uncertainties. 

With these promising results in mind, its usefulness in actual experimental practise 
needs to be demonstrated. Furthermore, one could widen the theoretical scope of our 
considerations towards higher-dimensional systems, as well as multipartite ones. 
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